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Quantum chaos and QCD at finite ctiemical potential 
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We investigate the distribution of the spacings of adjacent eigenvalues of the lattice Dirac operator. At zero 
chemical potential fi, the nearest-neighbor spacing distribution P{s) follows the Wigner surmise of random matrix 
theory both in the confinement and in the deconfinement phase. This is indicative of quantum chaos. At nonzero 
chemical potential, the eigenvalues of the Dirac operator become complex. We discuss how P{s) can be defined 
in the complex plane. Numerical results from an SU(3) simulation with staggered fermions are compared with 
predictions from non-hermitian random matrix theory, and agreement with the Ginibre ensemble is found for 
(1^0.7. 
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1. Introduction 

The spectrum of the QCD Dirac operator, both 
in the continuum and on the lattice, has several 
universal properties. By "universal" one means 
"independent of the details of the dynamics" , e.g., 
independent of the precise values of the simula- 
tion parameters on the lattice. Such universal fea- 
tures can be described by random matrix theory 
(RMT). In this contribution, we are concerned 
with the eigenvalue fluctuations in the bulk of 
the spectrum on the scale of the local mean level 
spacing, measured by the distribution P{s) of the 
spacings s of adjacent eigenvalues. We will con- 
sider gauge group SU(3) and staggered fermions 
which are related to the chiral unitary ensemble 
of RMT. At zero chemical potential ^, all Dirac 
eigenvalues are purely imaginary, and it has been 
shown in lattice simulations that P{s) agrees with 
the Wigner surmise of RMT, 
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both in the confinement and in the deconfine- 
ment phase 0]. This result implies that the 
Dirac eigenvalues are strongly correlated, and is 
indicative of quantum chaos, according to the 
conjecture by Bohigas, Giannoni, and Schmit 
j2|. In contrast, quantum systems whose classical 
analogs are integrable obey a Poisson distribu- 
tion, Pp(s) = e"". 

* Poster presented by R. PuUirsch 



Here, we focus on nonzero chemical potential. 
For /i 7^ 0, the lattice Dirac matrix generalizes to 
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with the Kawamoto-Smit phases 77. The eigen- 
values of this matrix are complex. In this case, 
P{s) represents the spacing distribution of near- 
est neighbors in the complex plane. For each 
eigenvalue zq one has to find the eigenvalue Zi for 
which s = \zi — zo\ is smallest, with a subsequent 
average over zq. This definition assumes that the 
spectral density is constant over a bounded region 
in the complex plane (and zero outside). Gen- 
erally, this is not the case so that an unfolding 
procedure must be applied, see Sec. |. However, 
the spectral density of the so-called Ginibre en- 
semble of RMT, where real and imaginary parts 
of the eigenvalues have the same average size, is 
constant inside a circle and zero outside, respec- 
tively Q. In this case, P{s) is given by [Q 

Pg{s) = epics) , (2) 

with 
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H = 0.l /i = 0.4 /i = 0.7 Ai = 1.0 ^ = 2.2 
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Figure 1. Complex eigenvalues of the Dirac operator at various values of /i for a typical equilibrium 
configuration of QCD with staggered quarks without chemical potential (horizontal axes = real parts, 
vertical axes = imaginary parts, in units of 1/a). 



In contrast, the Poisson distribution in the com- 
plex plane, representing uncorrelated eigenvalues, 
becomes 

Pp{s) = ^se-^^\ (3) 

In the following, we will study the Dirac spec- 
trum on the lattice at various values of /i 7^ and 
compare the resulting P(s) with Eqs. (^) and (|^). 

2. Analysis of complex spectra 

While the quenched approximation at /i 7^ is 
unphysical ^ , there is currently no feasible solu- 
tion to the problem of a complex weight function. 
Therefore, we generated the gauge field configu- 
rations at /J, = and added the chemical potential 
to the Dirac matrix afterwards. The simulations 
were done with gauge group SU(3) on a 6'^ x 4 lat- 
tice using /3 = 5.2 and Nf = 3 flavors of staggered 
fermions of mass ma = 0.1. For each parameter 
set, we sampled 50 independent configurations. 

The eigenvalue spectrum is shown in Fig. |l| 
for five different values of /i. The size of the 
real parts of the eigenvalues grows with /i as 
expected. While the spectrum is confined to a 
bounded region in the complex plane, the spec- 
tral density is certainly not constant. Therefore, 
the spectrum has to be unfolded. In one dimen- 
sion, unfolding is a local rescaling of the eigen- 
value density such that the density on the un- 
folded scale is equal to unity. To the best of 
our knowledge, unfolding in the entire complex 
plane has not been considered yet (in Ref. ||], 
unfolding was done in strips perpendicular to the 
real axis). We proceed as follows. The spectral 



density has an average and a fluctuating part, 
pix,y) = p^^{x,y) + pfi{x,y). Unfolding means 
to find a map z' — x' + iy' = u{x,y) + iv(x,y) 
such that piivix',y') = 1. Since the probability 
has to be invariant, pnv{x' ,y')dx'dy' — dx'dy' = 
P'Av{x,y)dxdy, we find that pa_^{x,y) is the Jaco- 
bian of the transformation from {x,y) to (x',y'), 
Pav(a;,y) = \dxudyV - dyud^^v] . Choosing y' = 
v[x,y) = y yields pav(a;,y) = \dxu\ and, thus, 
x' = u{x, y) = dtpa_v{t, y) . Essentially, this is 
a one-dimensional unfolding in strips parallel to 
the real axis. For a fixed bin in y, Pa-v{x, y) is de- 
termined by fitting p{x, y) to a low-order polyno- 
mial. P{s) is then constructed from the constant 
unfolded density as explained in Sec. |^, normal- 
ized such that ds s P{s) = 1. 

Several remarks are in order, (i) We have 
checked spectral ergodicity. If only parts of the 
spectral support are considered, the results for 
P{s) do not change, (ii) If the spectral density 
has "holes" (see Fig. |l| for ^ = 1.0 and 2.2), 
we split the spectral support into several convex 
pieces and unfold them separately. This is jus- 
tified by spectral ergodicity. (iii) Unfolding each 
spectrum separately and ensemble unfolding yield 
the same results for P{s). (iv) The results for 
P(s) are stable under variations of the degree of 
the fit polynomial and of the bin sizes in x and y. 

3. Results and discussion 

Our results for P{s) are presented in Fig. ^. 
For small values of /u, we probe the regime of 
"weak non-hermiticity" |^ where the typical size 
of the real part of the eigenvalues is of the order 
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Figure 2. Nearest-neighbor spacing distribution of the Dirac operator with complex eigenvalues for various 
values of /i. The histograms represent the lattice data. The short-dashed curve is the Ginibrc distribution 
of Eq. (|), the dotted curve in the left plot the Wigner distribution of Eq. (|l|) , and the long-dashed curve 
in the right plot the Poisson distribution of Eq. (||) , respectively. 



of the average separation of the imaginary parts, 
see Fig. |l| (our definitions differ from Ref. |^ by 
a factor of i so that it would be more appropriate 
to speak of "weak non-anti-hermiticity"). Some 
analytical results are known |^ in this regime, 
but the spacing distribution is only given as an 
expansion of P{s, zq) for small s, where zq is the 
location in the complex plane (i.e., no unfolding 
is performed). Therefore, we have no closed an- 
alytical result to compare with, and presumably 
it would be very difficult to derive one. However, 
something can be said about the small-s behav- 
ior of P{s). For the Ginibre ensemble, we have 
Pg(s) for s ^ 1. In the regime of weak 

non-hermiticity, it is known that P{s) starts with 
a power smaller than three ||^. This is conform 
with our results for fi ~ 0.1 in Fig. ^ where the 
Wigner surmise Pw(s) of Eq. (0) for the corre- 
sponding hermitian matrix with = is inserted. 

Increasing /i, we observe stronger level repul- 
sion leading to quite nice agreement of the data 
for fj, = 0.7 with the prediction from the Ginibre 
ensemble. For this value of /x, we are in the regime 
of "strong non-hermiticity" where the Ginibre en- 
semble applies. 

For values of /i larger than 0.7, the lattice data 
deviate substantially from the Ginibre distribu- 
tion. The main question is if our result for /i = 2.2 
in Fig. H should be interpreted as a transition 
to Poisson behavior, or if the agreement of the 
data at this particular value of fi with Pp{s) of 
Eq. (H) is purely accidental. We have to keep in 
mind that we neglected fi in the generation of the 
equilibrium gauge fields. Therefore, we are not 



entitled to make any definite statements on pos- 
sible connections between the finite density phase 
transition (expected at smaller values of /i) and 
the observed deviations from Ginibre behavior. 
Nonetheless, these deviations are interesting and 
deserve further study. 

In conclusion, we have investigated the nearest- 
neighbor spacing distribution of the lattice Dirac 
operator with finite chemical potential. The data 
agree with the RMT prediction for the Ginibre 
ensemble for /i = 0.7 and show deviations for 
smaller and larger values of fi. This is somewhat 
different from the situation at /i = where the 
data agree with the RMT results [Q. 
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